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A STUDY OF THE STRESS DISTRIBUTION 
AROUND ELLIPTICAL HOLES IN A FLATE PLATE 


SUBJECTED TO IN-PLANE LOADS 


by 


Demosthenes HOZOS 


ABSTRACT 


A theoretical solution of the title problem was obtained 
and a computer program was developed to examine various geometric 
arrangements of elliptical holes in flat plates with different 
Loading conditicns. 


Graphs of theoretical predictions are presented, discussed 


and compared with experimental results of other authors. 
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NOMENCLATURE 
Distance between the foci of the ellipses 
Major and minor semiaxes of the elliptical hole 
Elliptical coordinates of a point P(x,y) 


Constant value of § of any point on the elliptical 
hole 


Stress functions 


Normal stress intensity in the tangential 
direction to § = constant 


Normal stress intensity in the normal direction 
to = constant 
Shear stress intensity 


The usual notation for stress intensities in 


Cartesian coordinates 


Dimension of the plate in the X direction 


Dimension of the plate in the yi direction 


Uniformally distributed load applied on the 
WwW. 
edge a 


x= axis 





and parallel to xX 


Uniformally distributed load applied on the 


cieex y= wy and parallel to y axis 


semlaxes of the confocal ellipse oe Eo passing 


PCxX,y) = P(8>2) 


Gmax or Graax 
Ga 62 


Scr Ct) or 


throveh a point 


Stress concentration factor 





ft 


Stress concentration factor = 


Sc ros Cl) 








H 
pe 


b/a. 


a/b 


my oe 


when loading parallel to the minor axis of the 
elliptical hole. 


when loading parallel to the major axis of the 
elliptical hole. 








INTRODUCTION 


The behavior of thin flat elastic plates perforated with 
circular holes is well. known for in-plane loading. 

Theoretical methods and experimental results have been 
published by many authors. Refs. (1-7) 

The case of one elliptical hole in an infinitely large 
plate is discussed in reference [1] for several loading conditions 
including the case where the major axis of the ellipse is inclined 
to the axis of the axial load. 

Durelli and Murray [10] verformed photcelastic tests on thin 
flat plates loaded with various combinations of biaxial stresses. 

Durelli, Parks and Feng [9] examined cases of rectangular 
meewees Wien one e€llipticel hole under axial Leading, again Using 
the photoelastic method. 

The purpose of this study is to approach the problem 
thecretically and obtain results for tne case of finite thin 
rectangular plates perforated with one elliptical hole and 
subjected to uniaxial or biaxial loading. 

A computer program was developed and the published results 


of available experiments were used for comparison. 








THEORY 


The use of elliptical coordinates simplifies the mathematical 
analysis of the stresses around an elliptical hole in a flat plate 
loaded in the plane of the plate. 

The equations of transformation from Cartesian to elliptical 


coordinates are 


c-cosh 5. cos 7 (1) 


*« 
tl 


c-sinh §. sim 9 (2) 


~< 
i 


The conditions 3 = constant and 1 = constant are represented 
by a family of confocal ellipses and a family of confocal hyperbolae 


respectively, that are orthogonal and may be expressed as 


x? ae (3) 
‘Ccosh*€ : c* sink’€ 7 
2 
<< ae _ ——— = 4 (4) 
c* <os*1 Go Sin 7) 


Figure fa represents the elliptical coordinates 5,7 
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Using the notation 
i i (coshag - COSKT) ) OD 
R 


the stress components can be expressed in elliptical coordinates in 


terms of a stress function D as 


mm me me eS 


6y = 1 Db  sinh® cosh Pier 
J* as? ysintn ¢sinh*e) OF 


_ sin? cosy od (6) 
1° Gin'n ¢sinh=) 7 


Gr_ 1 DP __sinncosy 9 
Se o7n* J (sin’n +sinh §) ov 


sinh§S cosh& OD co 
LS, +sinh®) OF 





Go- L 2h, _smrecsy oP 
38 Onods J* (sin'n 4 sinh’s ) Ok 
sinh$ coshé oO (3) 


ae (sin’n +Sinh’s ) on 








See 


It is easy for one to verify that the expressions on page 
538 in [4] reduce to equations (6), (7) and (8) above. 

It might be noted here that the following assumptions have 
been made in the foregoing equations: 

i.Isotropic elastic material 

i, han fle place 

iii. Infinitesimal deformations 

The stress functions (da~-j) assumed here are the same as 


those which Coker and Filon used in reference [}] 


ec cos fom] 4 ial cos [(41)9] (92) 

xe = ies [on -10°9] i i cos[(n+1)7] (So ) 

x. J sinfcn-007] A BD em omtye] (9¢) 

.a= a sin| (n-1) 7] i lh sin] n+7| (9a) 

Xs - cos NY (9e) 

ae cosa (9£) 
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x7 = Ee sin oy (9g) 


a 
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Xa= - sin ny (9h ) 
Kqa= 5 (9i) 
io — uh (93) 


All the above functions satisfy the biharmonic equation for 


interger (or zero) values of E I. 


Equations (9a) throuch (9j) can be combined to give a general 


expression of the stress function 


ar. ke cngiys ; Cn-1)F 
Ga dot + Do +7. { Cn le COs CH-07) + 2 Cos (n4i)n| 4 
nN=90 


~(n-t ~(n-1)§ 
i dr fe a Os) (CH= ili} + 2 ° COS (R417 | 4 


ngs 


“hs en| € Sin (n-1)n + as sin cn4t)n | ‘ 
(nee. ~(a-)§ 
l: fn} e Sin (n-1)y +e sin ntiy | i 


n§ ~ne 
+ Se cosny 4 hye cosnn 4 


é ne ns . 
+ jn @ Simny tne sin ay (10) 





= eee 


When examinig the series with n=O and n=l then for 
Mimicte stresses and displacements 16 is evident tnacetire 
@octficients Dot Mos Igo K aot, must equal zero. 

The stress function given by equation (10) can be 
substituted in equations (6),(7),and (8) for the stresses. 

However am the present analysis,only Word tre ic oncom. 
that maintain symmetry about both the x and y axes are 
considered anc therefore the coefficients of the terms that 
Ppmeoeauce nonsymmetric functions of iy, for different Vvanieamoman 
were set equal to zero. 


Finally the stresses from (6),(7) and (8) have the form: 
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or Ge 


However equations (11), (12) and (13) are rearranged, in Appendix 1, 
in a more suitable form for translation into computer language. 
The normal and shear stress components in Cartesian coordinates 


can be computed from the equations, 


Gx - Sy oe ee f Ge sinc ~ Tn SIN 2x (14) 
Care - (G: -Gn) SIN RX ~ Ten Cos RK (16) 
where cos Va_cosh§ sing 


SON Gis se 


sina — V2_sinh¥ cos (18) 


J cosh es -cos en 


Equations (14), (15) and (16) are derived as usual from 


equilibrium conditions on an element .(rig. 1b) 
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BOUNDARY CONDITIONS 


Equations (11), (12) and (13) satisfy the required conditions 
of symmetry and antisymmetry refered to in the previous section. 
In addition to this Ge and Csm must vanish at &=35o. Page 198,| 4 | 


Therefore at each point on the elliptical hole we must satisfy 


Sr = © Be Ge (19a) 
Ty7 =O 8305 F=Fo (19 ) 
For the case of one elliptical hole in an elastic finite 


flat rectangular plate loaded biaxially we must satisfy the following 


two conditions (20) and (21), in addition to (19) 


Txy 220) (20a) 
Gy ere a bf = haa (206) 
A 
exe O (Ala) 
Gx = Gi x = Wx (21 b ) 
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Due to the symmetry about both the x and y axes it is 
enough to examine only one quadrant of the plate. 

Figure 2 illustrates all the cases under consideration. 
The numbers appearing in parentheses are the numbers of the 
equations of the boundary conditions applied for each particular 


case. 
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NUMERICATL, METHOD 


Initially the number of unknown constants appearing in equations 
(11) up to (13) was determined by an arbitrary truncation of the series. 

Later the behavior of the series was examined for several cases 
and it was realized that the series doesn't always behave well. Similar 
convergence problems for functions involving elliptical coordinates were 
pointed out in reference [6] 

In this case, the reason is believed to be the very rapid convergence 
of the coefficients of Dn and Hn after the fifth or sixth term, resulting 
in avery rapid divergence of the coefficients of Cn and Gn for all of 
the cases that were examined. 

Finally twenty terms plus the constant term of the series were 
e7osen. 

Using the Overdetermined Collocation method (Ref.8 ) the following 
eases were examined: 

ae- Six fixed points on the inner boundary and another six fixed 

points on the outer boundary were selected to satisfy the 
boundary conditions. (Fig.2d) 
b.e~ Decreasing by one the number of points on the ellipse with 
no change on the outer edges produced no difference in the 
numerical results obtained in case (a). 

C.~ With six points on the ellipse and eight on the outer 
boundary, the shape of the curve for the tangential stress 
distribution around the elliptical hole was changed to a 


flatter one and it was realized that the boundary conditions 





ene 


were not satisfied as well as before due to the least 
squares process involved in the numerical method. 

d.- The case of five points on the inner boundary and eight 
on the outer one was examined and it was found that 
the error introduced was smaller than in case (c), but 
not as small as in cases (a) and (b). (Verifying the 
satisfaction of the boundary conditions) 

It might be noted here that a smaller number of points than 

in case (a) was used by Christiansen [ 15 J] for a similar problem. 

For the case of axial loading parallel to the minor axis 

of the ellipse and for values of 4% <0.5 (wide plate) the 
expression & approaches infinity and no results could be obtained 
(Appendix 1). 

The same effect happens for the case when plates with dAy>05 

are loaded parallel to the major axis of the elliptical hole. 

From the boundary conditions a set of M linear algebraic 

equations in WN unknowns is formed, where MON . 


In matrix notation 


AX 


i 
es) 


(22) 


where A is the M by N matrix of the coefficients of 


the unknown constants X. 


By means of a Least-squares method the normal equations 
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are obtained and the system in matrix notation has the form 


AAX = AB (23) 


These normal equations are solved by elimination using 
largest pivotal divisor [Ref. 10 ] 

A computer program was developed and written in FORTRAN IV 
to execute all the required steps and calculate the tangential 
stress Gy around the ellipse and the normal stresses along the 
axes of symmetry of the plate. 

Subroutines for matrix operation are taken from reference [ 12] 

Ail the necessary calculations for the problems under 
consideration are carried out by the same computer program by 
transfering control to the proper boundary conditions. 

The required input information for each plate to be 
examined, must be typed on one card and many plates can be 
analyzed the same time. 

Details for the input cards are given on comment cards in 
the program. 

Additional details about the problem formulation are given 


mo Appendix lL. 
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RESULTS 


Figures 3 up to 31 are typical plots of results obtained 
from the computer program. Whenever possible, experimental data 
are indicated. Most of the figures are selfexplanatory but additional 
information, if necessary, is given in the discussion of each figure. 

Figures 3 up to 14 show the stress distribution around the 
poundary of the elliptical hole in finite plates for several values 
of He and. q° 

Tt should be noted that in cases of loading parallel to the 
minor axis, the tensile stresses predicted are 2% to 17% higher 
than experiments reported in [ 9] while the much smaller compressive 
stresses are 40% to 604 lower. 

When the external. loading is parallel tc the major axis of the 
ellipse the agreement with [9] lies between 5% and 25%. 

The above differences in compressive stresses, when Loading 
parallel to the major axis are of minor importance from an engineering 
point of view because, these stresses are very small compared to the 
maximum tensile ones. 

Figures 15 up to22 show the maximum tensile stresses (SCF and 
Scr!) as function of q for several values of WM . (Finite plate 
with one elliptical hole). 

In figure 25 several cases of biaxial loading are plotted. 


Tt was not possible to examine the same olate used in [10] because 


one 7 limitations discussed previously. Other plates were examined 





aA 


and the'definitions from [10] were used. 

The curves appear to have the same shape. 

Figures 2h up to 31 show the 67 and Gy distribution around the 
boundary of the elliptical hole and the x axis respectively for typical 


uniaxial and biaxial loading. 
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CONCLUSIONS 


The stress functions assumed have been shown to give results 
that are comparable with available experimental data and therefore 
some confidence exists in the theoretical predictions for cases 
where no experimental result3 exist. 

The stress concentration factors for plates perforated with 
elliptical holes under in-plane loading, appear to be much higher 
in comparison with Similar plates perforated with circular holes. 

Certain disagreements have been noticed between measured 
and predicted stresses in the low level compressive stress range, 
while the predictions for the high level tensile stresses are 


somewnat higher than the measured ones but still in close agreement. 
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RECOMMENDATIONS 


The following topics are recommended for further investigation: 

a.- Postulation of similar stress functions in terms of 
hyperbolic functions instead of positive exponentials, in 
association with other functions decaying rapidly as 
increases in order to be more suitable for computer use. 
Timpe in [13] suggests such stress functions. 

be= Study of the problem for cases of loading that produce 
antisymnetrical conditions (shear forces applied on the 


edges of the plate for example), will be rather easy to handle. 
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APPENDIX 1 


Equations (11), (12) and (13) may be rearranged in the form of (2), 
(25) and (26), respectively. Also equations (14), (15) and (16) are 


rearranged in the form of (27), (28) and (29) respectively. 
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In the above expressions for ©x, Sy , Txy 
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It is obvious that as )+0 (i.e. XKWk) dD'» a!’ and therefore 
SS 

Analyzing the case of figure ed, for example, the use of 
equations (24) and (26) according to (19a) and (19b) generates a set 
of le equations, if six fixed points are selected on the boundary of the 
(quadrant) ellipse 

Another set of six equations is generated by using (27) and (29) 
according to (21), if three points are selected on x=3, 

Finally six more equations are generated by using (28) and (29) 
according to (20), if another three points are selected on yea3e 

It should be noted that the number of unknown constants was 
selected to be 21. 

The system obtained this way has 24 equations with 21 unknowns 
and it is solved by following tne vrocedure described in section 
NUMERICAL METHOD. 

The matrix A! is formed and A and B are premultiplied by 


AT and from the solution of the system 


ATAX = AT B 


we obtain the values of the unknown constants. 
Summarizing the above procedure, we have chosen 2 number of 
points greater than the number of unknown constants and we required 


that the boundary conditions be satisfied as nearly as possible at 





ar eae 


these points (Least square sense) Ref. [ 14]. 

Because of the nature of the assumed series solution, the number 
of unknown constants is always odd when truncating the series at any 
value of n. Therefore the system of cquations is always overdetermined 


at least by one equation. 
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